Shake-up and shake-ofF excitations with associated 
electron losses in X-ray studies of proteins 



Petter Persson ^ , Sten Lunell ^, Abraham Szoke 
Beata Ziaja ^'^'^ and Janos Hajdu ^ Q 

^ Department of Quantum Chemistry, Uppsala University, 
Box 518, S-751 20 Uppsala, Sweden 

^ Department of Biochemistry, Biomedical Centre, 
Box 576, Uppsala University, S-751 23 Uppsala, Sweden 

^ Lawrence Livermore National Laboratory, Livermore, CA 94551, USA 

^ Department of Theoretical Physics, Institute of Nuclear Physics, Radzikowskiego 

152, 31-342 Cracow, Poland 

^ High Energy Physics, Uppsala University, P.O. Box 535, S-75121 Uppsala, Sweden 
Corresponding author: 

Janos Hajdu, Department of Biochemistry, Biomedical Centre, 

Box 576, Uppsala University, S- 75123 Uppsala, Sweden 

Tel:+4618 4714999, Fax:+4618 511755, E-mail:liajdu@xray.bmc.uu.se 

Manuscript information: 

Number of pages (double-sided version with tables and figure legend on separate pages): 
19 , Number of figures: 1, Number of tables: 2 

Submission information: 

Diskette of the manuscript included, manuscript written in LaTeX, ASCII file ex- 
tracted; hard copy double-spaced 



-mailipetter. persson@kvac.uu.se, sten.lunell@kvac.uu.se, szokel@llnl.gov, 
ziaja@tsl.uu.se, hajdu@xray.bmc.uu.se 



Abstract 



Photoionisation of an atom by X-rays usually removes an inner-shell electron 
from the atom, leaving behind a perturbed "hollow ion" whose relaxation may 
take different routes. In light dements, emission of an Auger electron is com- 
mon. However, the energy and the total number of electrons released from the 
atom may be modulated by shake-up and shake-off effects. When the inner-shell 
electron leaves, the outer-shell electrons may find themselves in a state that is 
not an eigen-state of the atom in its surroundings. The resulting collective exci- 
tation is called shake-up. If this process also involves the release of low energy 
electrons from the outer shell, then the process is called shake-off. It is not clear 
how significant shake-up and shake-off contributions are to the overall ionisation 
of biological materials like proteins. In particular, the interaction between the 
out-going electron and the remaining system depends on the chemical environ- 
ment of the atom, which can be studied by quantum chemical methods. Here we 
present calculations on model compounds to represent the most common chem- 
ical environments in proteins. The results show that the shake-up and shake-off 
processes affect about 20% of all emissions from nitrogen, 30% from carbon, 40% 
from oxygen, and 23% from sulphur. Triple and higher ionisations are rare for 
carbon, nitrogen and oxygen, but are frequent for sulphur. The findings are 
relevant to the design of biological experiments at emerging X-ray free-electron 
lasers. 

Keywords: X-rays, photoionisation, shake-up, shake-off. Auger emission, radiation 
damage, peptides, proteins 
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Computer simulations (Neutze et ai, 2000) show that ultrashort and high-intensity 
X-ray pulses, as those expected from presently developed free-electron lasers (Winick, 
1995; Wiik, 1997), may provide structural information from large protein molecules and 
assemblies before radiation damage destroys them. Estimation of radiation damage as a 
function of X-ray photon energy, pulse length, integrated pulse intensity and sample size 
was obtained in the framework of a radiation damage model, where the effects of atom- 
photon and ion-ion interaction were taken into account. Photons of 1 A wavelength, 
corresponding to a photon energy of about 12 keV, interact with atoms mainly via the 
photoelectric effect (Dyson, 1973) (for carbon the photoelectric cross-section is ~ 10 
times higher than the corresponding elastic cross-section at this wavelength (Hubbel 
et al, 1980)), and thus the photoelectric effect is the main source of radiation damage 
with X-rays. 

Photoionisation may proceed either through the ejection of an outer-shell electron 
or through the ejection of an inner-shell electron (Dyson, 1973). Outer-shell photo- 
events eject a single electron from the atom with an energy equivalent to the energy of 
the incoming photon minus the shell binding energy and the recoil energy. With photon 
energies of around 12 keV (or about 1 A wavelength), outer-shell events are rare in 
carbon, nitrogen, oxygen and sulphur, and represent less than 5% of all photoionisation 
events. 

A more frequent type of photoionisation with X-rays involves the ejection of an 
inner-shell electron from the atom (around 95-97% of photoionisations remove a K- 
shell electron from carbon, nitrogen, oxygen and sulphur). Typical K-hole life times 
are of the order 1 — 10 fs (Krause & Oliver, 1979) , and the hollow ion relaxes through an 
electron falling from a higher shell into the vacant hole. In light elements, the energy 
of this electron is given to another electron which is then also ejected from the atom 
through the Auger effect. Core ionisation, however, constitutes a strong perturbation 
of the molecule, and may be accompanied by significant electronic effects (Siegbahn 
et ai, 1969). Firstly, the valence electrons often relax significantly to compensate 
for the presence of the positive core hole. The departing photoelectron can interact 
with these relaxing electrons, and lose kinetic energy in the process, thus forcing the 
system into an excited state (a process called shake-up). In some cases, the excitation 
may result in the ejection of low energy outer-shell electrons from the atom (a process 
called shake-ojf). The multiple excitation lifetimes for light elements arc comparable to 
core-hole lifetimes. The interaction between the out-going electron and the remaining 
system depends on the chemical environment of the atom, and a description of these 
interactions requires quantum chemical calculations. The following key processes need 
to be taken into account: 
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(i) PHOTOEMISSION FOLLOWED BY SINGLE AUGER DEGAY. 

Ejected K-shell electrons with the highest possible kinetic energy correspond to a situa- 
tion where the remaining system is left in the ground state, i.e. the difference between 
the incoming photon energy and the energy of the ejected photoelectron equals the 
chemical binding energy and the recoil energy. In light elements like carbon, nitrogen, 
oxygen and sulphur, such clean photoemissions are followed by Auger emission, and 
the atom becomes doubly ionised (Krause & Oliver, 1979). In heavier elements X-ray 
fluorescence dominates. 

(ii) SHAKE-RELATED PROGESSES. When photoionisation ejects an inner-shell 
electron so fast that the outer-shell electrons have no time to relax, the situation gets 
similar to beta decay, whereby the nuclear charge suddenly increases by one unit. 
Quantum mechanics describes such a state as a superposition of proper eigen-states, 
that include states where one or more of the electrons may be unbound. Interactions 
between the departing photoelectron and the electrons left behind may reduce the 
kinetic energy of the photoelectron, and deposit energy into the system. The pertur- 
bation is called shake-up if they refer to an excitation in the final system, or shake-off 
if the result is the loss of one or more outer shell electrons from the ion. The relative 
contributions from these two processes have been found to be comparable in noble 
gases (Svensson et al, 1988; Armen et al, 1985; Wark et al, 1991). The following 
shake-related phenomena need to be considered: 

(a) Photoemission accompanied by shake-up excitation with Auger decay. This 
process reduces the energy of the photoelectron shghtly (about 10-40 eV), and the 
energy difference can either be absorbed by the atom or added to the energy of the 
Auger electron. At the end, the atom becomes doubly ionised. 

(b) Photoemission accompanied by a shake-off event. In the vicinity of a hole, a 
vacancy and a free electron are created. The shake electron has around 10-100 eV 
energy, and the atom becomes doubly ionised. 

(c) Photoemission accompanied by double Auger decay, which may proceed via 
different mechanisms, and may result in the triple ionisation of the atom. These mech- 
anisms include Auger cascading, single Auger decay combined with shake-off emission, 
and virtual inelastic scattering (for a full description see Amusia et al, 1992). In this 
case, the energies of the two ejected electrons are asymmetrically distributed. For in- 
stance, for Ne (transition Is^^ 2s^^2p~^ + ?i + 92) the total available kinetic energy 
is about 650 eV, and the most probable case is shaking off a slow electron and Auger 
emission of a fast electron with Eshake-off ^ E Auger (Amusia et al, 1992). 

In the initial damage model of proteins (Neutze et al., 2000) only outer-shell pho- 
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toionisation and inner-shell photoionisation with a subsequent Auger emission were 
taken into account. In the present paper we focus on photoionisation mechanisms, 
including shake-up with a subsequent satellite photoionisation (case a) and shake-off 
(cases b and c), which under certain conditions may produce multiple ionisations in 
the atom (case c). We investigate in detail what effects these processes may have on 
model compounds in order to assess radiation induced damage processes in biological 
samples. We calculate the different shake contributions for carbon, nitrogen, and oxy- 
gen atoms in a model peptide (see Fig. la). Shake effects for sulphur are estimated in 
separate calculations for the three most common chemical environments of sulphur in 
proteins (Cys, Met, Cys-S-S-Cys; Figs. Ib-d). 

Results 

Fig. 1 shows the model compounds used in the calculations. In a recent polymer study 
(Nakayama et al, 1999), experimentally observed differences in intensity for ester and 
carbonyl oxygen main lines could be explained directly from consideration of the main 
line intensities, where the main line intensity was estimated from the overlap obtained 
from the n = term in equation (|l]). This has the clear advantage for inherently large 
molecules, such as polymers or polypeptides, that main-line losses can be estimated 
without consideration of all the shake states to which intensity is distributed. For large 
molecules, the number of shake states goes beyond that which can readily be calculated 
at the configuration interaction level of theory. In a Gly-Gly -Gly tripeptide (Fig. la), 
the central glycine unit can be expected to display a shake spectrum which is similar to 
that of a residue within a longer polypeptide chain. To investigate effects arising from 
the truncations of the polypeptide chain, shake effects were calculated for all atoms 
in the model, including the atoms from the terminal carboxyl and amino functional 
groups. 

There are six carbon atoms in the peptide model, showing very similar overlaps 
between the initial unionized state, and the core ionized ground state (cf. Table 1). 
This overlap is ~ 0.85 for all carbons, except for the carbon at the terminal carboxylic 
acid group, which has an overlap of ~ 0.87. The intensity loss from the main line is 
therefore close to ca. 30% for all carbon atoms in these models (cf. Table 1). 

The oxygen shake contribution is generally the largest. In particular, the carbonyl 
oxygens display large shake effects, in agreement with the general behaviour for organic 
molecules. Here, strong shake is found for the two carbonyl oxygens in the peptide 
links, as well as for the double-bonded oxygen in the terminal carboxylate functional 
group. The overlap for these three oxygens is ~ 0.78 (cf. Table 1), corresponding to a 
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total shake intensity loss of ca. 40% (cf. Table 1). The terminal hydroxyl oxygen has 
a slightly larger overlap of ~ 0.82, corresponding to the main line intensity loss of ca. 
33%. 

There are two nitrogens in the peptide links, and one in a terminal amino group. 
The shake effect for all these atoms is similar: the overlap is ~ 0.90 (Table 1) and the 
main line loss is ca. 20% (cf. Table 1). Nitrogen shows the smallest shake contribution, 
and the results agree well with data on a polyimide polymer (Nakayama et al., 1999), 
which also showed little nitrogen shake effect. 

Sulphur shake effects have been calculated for three different common chemical 
environments (Fig. Ib-d and Table 2). The largest loss was found for a disulphide 
bridge (R-S-S-R, 25 %), and the smallest one for the thiol terminal group (R-SH, 21 
%). 

The initial models used for the sulphur calculations were smaller than the models 
used for the other atom types. This could, in principle, be a cause of concern, as 
these type of calculations are known to underestimate shake contributions in extended 
systems. In order to test the stability of the numerical results, additional calculations 
were performed on models with increasing chain lengths ( CH^ — SH, CHs — CH2 — SH, 
and CHs — CH2 — CH2 — SH ). The results indicated negligible drift (not shown). 

Discussion 

Radiation damage prevents the structural determination of single biomolecules and 
other non-repetitive structures (like cells) at high resolutions in classical electron or 
X-ray scattering experiments (Henderson, 1990; Henderson, 1995). Analysis of time- 
dependent components in damage formation suggests that the conventional damage 
barrier can be substantially extended at extreme dose rates and ultrashort exposure 
times (Neutze et al, 2000; Hajdu, 2000; Hajdu et al, 2000; Hajdu & Weckert, 2001; 
Ziaja et al, 2001). A quantitative description of damage formation and a detailed 
analysis of the ionisation dynamics of the sample are crucial for planning experiments 
at future X-ray lasers. Results described in this paper give the first assessments of 
shake contributions to the ionisation of carbon, nitrogen, oxygen and sulphur atoms 
within the most common chemical environments in a protein molecule. In previous 
applications, the calculated shakeup intensities were usually a few percentage units 
too high, so we judge the present results to provide upper bounds to the total shake 
intensities. 

Data for noble gases show that shake-up and shake-off effects were of comparable 
magnitude (Svensson et al, 1988; Armen et al, 1985; Wark et al, 1991). Calculations 
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on carbon, oxygen and nitrogen by Mohammedcin et al., 1993 and El-shemi & Hassan, 
1997 show that photoemission accompanied by double Auger decay has nearly zero 
probability for these light elements. Such triple ionisations can thus be neglected in 
the ionisation dynamics of C, O, N compounds. Sulphur, on the other hand, behaves 
differently, and the probability of multiple ionisation (n > 2) after the ejection of a 
K-shell photoelectron is larger. The average charge left in sulphur after a single K-shell 
photoionisation is estimated to be about 4 (Mohammedein et al., 1993). For instance, 
for S: P(2) ~ 5%, P(4) ~ 34%, P(6) ~ 3% (Mohammedein et al, 1993; El-shemi & 
Hassan, 1997). This implies that for sulphur, multiple ionisation mechanisms should be 
taken into consideration in order to obtain more accurate predictions for the ionisation 
dynamics of proteins. 

The energy of the photoelectrons, which emerge from shake-up processes is similar 
to the energy of photoelectrons released during outer-shell ionisation (Nakayama et al, 
1999). In these calculations the deviation was within 10-40 eV. As a consequence, 
these electrons can leave a submicroscopic sample in the same manner as outer-shell 
photoelectrons. The inelastic mean free path of these electrons is of the order of a few 
hundred Angstroms (Ashley, 1990 and references therein; Ziaja et al, 2001), and thus 
the damage by the departing electrons will only have to be considered in large samples 
where such electrons may become trapped and deposit energy. 

Low energy shake-off electrons will behave similarly to Auger electrons, and may 
cause secondary ionisation in the sample. A detailed analysis of secondary electron 
cascades elicited by low energy electrons can be found in Ziaja et al, 2001. 

Materials and methods 

The average total shake contributions were estimated from quantum chemical (QC) cal- 
culations, using the semi-empirical INDO/S-CI program ZINDO, developed by Zerner 
and co-workers (Ridley & Zerner, 1976; Bacon & Zerner, 1979; Zerner et al, 1980), 
with standard parameterisation. A single determinant description was used for the 
initial (unionized) state of the system, while the excited core-ionized states were ob- 
tained from Configuration- Interaction calculations with Single excitations (CIS). The 
equivalent core approximation was used to represent the core ionized species. 

The calculations are based on the sudden approximation (Aoberg, 1967) valid for 
high energy photons, in which the intensity of a given shake line is proportional to the 
square of the overlap (^ol^n) '■ 



Pn = |(*o|$n)|', 



(1) 
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where is the wave function of the N-1 remaining electrons in the neutral molecule, 
and is the (N-l)-electron wave function of the nth state of the ionized system. 
These intensities were calculated with the SHAKEINT (Lunell, 1987; Lunell & Keane, 
1988) program package. This method has been successfully applied to a wide range of 
systems, including organic molecules (c.f. Lunell et al., 1978; Sjogren et al, 1992), C60 
(Enkvist et al., 1995), polymers (PMDA-ODA polyimide) (Nakayama et al., 1999), and 
adsorbates on metals (CO adsorbed on Cu(lOO) surfaces) (Persson et al, 2000). 
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Table 1. Overlaps ((^o|*^o)) '^^'^ average total shake electron intensity (1 — |(^o|*^o)P)- 
Values were calculated for all carbon, oxygen, and nitrogen atoms ocurring in the model 
polypeptide (Fig. la). 



Atom 


Overlap 


Average shake intensity 


CI 


0.853 


27% 


C2 


0.858 




C3 


0.850 




C4 


0.857 




C5 


0.853 




C6 


0.869 




01 


0.776 


38% 


02 


0.773 




03 


0.785 




04 


0.823 




Nl 


0.902 


19% 


N2 


0.896 




N3 


0.897 





Table 2. Total shake contributions (1 — |(^'o|$o)P) calculated for sulphur in three dif- 
ferent chemical environments. 



Molecule 


Shake intensity 


C H^C H2 — SH 
C H^C H2 — S — C H2C 
C H^C H2 — S — S — C H2C 


21% 
22% 
25% 


Average 


23% 



Figure 1. The model polypeptide for calculating the shake contributions for carbon, 
nitrogen, and oxygen atoms (Fig. la). Shake effects for sulphur were estimated in 
separate calculations for the three most common chemical environments of sulphur in 
proteins (Cys, Met, Cys-S-S-Cys; Figs. Ib-d). 
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